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^ ■ Abstract 
O ■ 

^ . We study the phase structure of a 4D complex scalar field theory with a 

Q \ potential V{^) = I ~ P a^t zero and at finite temperature. The model is 

analyzed by mean-field and Monte Carlo methods. At zero temperature the 
theory falls in the universality class of the 4D Ising model when varying A. 
The situation is less clear-cut for variations w.r.t. A at large temperatures and 
variations w.r.t. temperature at a given value of A. We observe temperature 
p) • independence of the mass of the first excitation. 

D ■ 

1 Introduction 

(N 
> 

^ ■ The most striking success of the theory of critical phenomena is its ability to explain 

OO . the coincidence of critical exponents in physical systems that, at first sight, seem to 

^ ! be vastly different 0. The behavior in the vicinity of continuous phase transitions 

\ can be explained in terms of just a few fundamental properties of the system under 

consideration: The dimension of spacetime and the symmetries of the associated 
Lagrangian. As a consequence, the tremendous variety of physical systems exhibiting 
critical behavior is classified by this universality. 

It thus is of a certain interest to study the nature of phase transitions in a system 
^ \ whose Lagrangian possesses unusual symmetries. Consider the field theory defined 

\ on four dimensional Euclidean spacetime by the Lagrangian density 

£($,9.$) = l{(a,$*)(a,$) + V(<l>)} (1) 
where the potential for the complex scalar field $ is given as 



o 



A<l> 



2 



. 2 1^,2 , _ g ,4 (Re$f-(Im$)^ 



^Confer e. g. [Tj for a review on this subject. 
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with A being a constant of dimension mass. This model has been proposed in 
the framework of a more general effective theory as described in [21 [3], H] and we 
refer the reader to these papers for details. A contour plot of the potential (E]) is 
shown in Fig. [TJ Notice that V{^) is symmetric under both of the following discrete 
transformations: $ t-^ — $ and $ i-^ $* or, equivalently, arg($) arg($) + tt 
and arg($) — arg($). Notice also that the shape of the potential suggests the 
existence of phase transitions. 

In this work we aim at obtaining the following information about the theory defined 
by the Lagrangian ([1]): First, we are interested in resolving the phase structure of the 
model in the parameter A and the temperature T. The latter will be introduced by 
an anisotropic coupling 7 in a lattice model. That is, we would like to identify the line 
in the A-7 plane which separates the phase with dynamically broken Z2- symmetry 
from the phase where this symmetry is respected by the ground state. Second, we 
intend to investigate the order of the phase transition and - if a continuous phase 
transition occurs - to determine its critical exponents. Finally, we also determine 
the mass of the first excitation at both zero and finite temperature. This is used to 
relate the bare asymmetry parameter 7 to temperature. 

The paper is organized as follows. In section [2] we define the lattice version of the 
continuum field theory given by the Lagrangian density ([H) and discuss the limiting 
cases A — > and A — 00. Subsequently, we carry out a mean-field analysis for the 
zero temperature case in section El which allows us to gain first insights into the 
phase structure of the theory. Section HI is concerned with a Monte Carlo study of 
the model, and it constitutes the central part of this work. We first design and test 
suitable Monte Carlo algorithms in section 14.11 In the remainder of section HI we 
obtain the phase diagram, determine the masses and the physical anisotropy, and we 
carry out a finite-size scaling analysis to determine the critical exponents. Finally, 
we summarize our results in section [5l 



2 Lattice version of the model 

We discretize the field theory defined by ([T]) on a four dimensional lattice in Eu- 
clidean spacetime. The field at site n is referred to as where n is a lattice point. 
The lattice spacings along the temporal and spatial axes are denoted by and ag, 
respectively. We label unit vectors by a hat, e. g. f , j, and their length is one in 
units of the respective lattice spacing. In particular, n + j is the nearest neighbor of 
site n along the (spatial) j— axis and similarly for n + f. Following common con- 
ventions, Roman indices refer to spatial coordinates while Greek indices run over all 
coordinates of the compactified Euclidean spacetime including the temporal axis. 
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Figure 1: A contour plot of the potential Eq. for A = 1.0. The axes are scaled 
in units of A. Note the minima at $ = ±A on the real axis, the saddle points at 
$ = ±z A and the pole at $ = 0. 



The discretized action then reads 



Here we have introduced an asymmetry parameter 7 (also referred to as bare 
anisotropy) between temporal and spatial couplings. This allows us to study the 
model at finite temperature along similar lines as in [HI |6] . We frequently decompose 
the field variables into their real and imaginary parts according to = + iUn- 
It will turn out to be convenient to define the following quantities: 



(3) 




's 



Sn := + $„-f ) + E i^nV, + ^n-^ 




(4) 
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At some points it will be convenient to write as $„ = RnC^^". In terms of these 
variables the lattice action reads 

r r 

S" = — ^ I -Rn (cos Qn+f COS 6^ + sin Qn+i sin e„) 

e„+-. cos e„ + sin e„+. sin e J (5) 

+ (3 + 7^ + + - a^A^ cos(2e.) } . 

Also, we adopt the following notation: For each dimensionful physical quantity (as 
for instance $, A etc.) we refer to the corresponding dimensionless quantity as mea- 
sured in appropriate powers of the spatial lattice spacing by the same symbol 
endowed with a hat (e .g. A etc.). 

Before addressing the dynamics of the model we analyze its behavior in the limits 
of weak and strong coupling A at zero temperature (7 = 1). For A — > we find 

^ = E(4 + Y)^n - E ^iRnSl-Sn + 0{A'), (6) 

n <l,n> 

where we have defined S„ := ( cos 6„ , sin 6„ ). The sum with the subscript < l,n > 
means summation over all pairs of nearest neighbors. This expression resembles the 
action of a free theory with a mass term. 

In the strong coupling limit A —> 00 all field configurations which do not have the 
field sitting in one of the two minima of the potential at ±A on the real axis are 
strongly suppressed. The remaining degree of freedom is the sign of the real part of 
the field, that is s„ := sign(Re$n). Its dynamics is governed by lowest-order terms 
in A in the Euclidean action which read 

5'a^oo:=-A^ SlSn, (7) 

</,n> 

Eq. (jTj) is nothing but the familiar Ising model in four dimensions. 



3 Mean-field analysis 

In this section we study the model at hand in the mean-field approximation (MFA). 
We restrict our analysis to the isotropic case 7 = 1. We find that the relevant 
symmetry is the Z2-symmetry relating the real minima of the potential. A phase 
transition occurs at a critical value Ac of the lattice coupling constant above which 
this symmetry is dynamically broken. In the following we shall stick to the "old 
fashioned" MFA (see, e.g. [7j) which is more intuitive. 
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Let us consider the contribution to the overall action of the fixed field variable $„. 
From Eq. (JSj) one finds that for 7 = 1 this contribution is given as 

S{^n} = -Rn ^ [cos 6n (^-Rn+A COS 6„+/i + Rn-f, COS Qn-fi 

+ sin e„, {kn+ji sin Qn+fi + -Rn-/i sin Qn-f}j 



+ (d+y)i?^ + ^-A^cos(2e. 



(8) 



= :V{kn) 

where d denotes the dimension. If not stated otherwise we set d = 4. 
The mean-field approximation amounts to decoupling each field variable $„ from its 
neighbors by replacing each of the 2 ■ d adjacent fields $n±/i by a mean-field me*'^, 
thus (dropping the indices 'n') Eq. ([H]) becomes: 

S{^) = -2rfm^ (cose cos + sine sin cu) -A^cos(2e) + 

= -2dmR cos(e - lu) - cos(2e) + V{R) . ^^"^ 

As a result the partition function factorizes and we only have to deal with a single 
integral: 

Z(A,m) = j d'<le-^(*) 

= j dRj c/e/?exp|+2c/mi?cos(e-cj)+A^cos(2e)- V(i?)} . 

(10) 

Selfconsistency then requires the expectation value of $ to equal the mean field, 
hence: 

-| roc /»27r 

me^^ = (<l) = — . / RdR dB e~^^^^®^ Re'^ 

Z{A,m) Jo Jo (11) 

=: (p{A,m) . 

This mean-field equation always has the trivial solution m = ( $ ) = 0. It can be 
shown from the symmetries of the integral that non-trivial solutions require 07 = 
(or vr) and u = ^ (or — f). For small A only the trivial solution m = exists. 
With increasing A a non-trivial solution appears at some value Ac, above which $ 
develops a non-zero expectation value m(A) ^ 0. Since this solution has a smaller 
free energy than the trivial on^ it replaces the latter, triggering a phase transition. 



^The trivial solution has in fact infinite free energy, but this is an artefact of the MFA for a 
potential singular at 0. 
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Figure 2: This figure illustrates the determination of A^; from the criterion 



0, the inset shows the vicinity of the zero in more detail 



0, 



solid curve). The dashed curve is the result for a; = |, showing a later onset. 



Some detail is given in appendix A.l. It turns out that the solution with i-i^ = f sets 
in at a larger A than the one with uj = and has a larger free energy than the latter, 
therefore it does not take over. This is physically sensible, since it corresponds to a 
mean field sitting on the saddle points of the potential and not at the minima. In 
the following we shall therefore take u = 0. See also Fig. [21 
We define 

^n(A,m):=^(A,m) (12) 
which can be obtained analytically (see appendix A.l). Then 

^(A.™) = i-^l^. (13) 

The onset of the symmetry breaking solutions is given by: 

^(Ac,0) = l. (14) 

This is displayed in Fig. [21 For later comparison to the Monte Carlo results we quote 
the numerical value 

AcMFA ^ 0.2920 . (15) 

It is straightforward to obtain further quantities, for example the critical exponent 
/3 which governs the behaviour of the field expectation value near the threshold: 

/A - Ar\(^ 

for A > Ac. (16) 
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Figure 3: Determination of the critical exponent (3 from the MFA. The plot shows 
\ d(inA) plotted versus In A. 



For this we expand ffT^ into a Taylor series in m. Making use of relation ffT2|) and 
noticing that -Fn(A, 0) = for odd n, we arrive at 

1 mF2(A,0) + im3F4(A,0) + O(m5) 

TTl — (li) 

2rf Fo(A,0) + ^F2(A,0) + 0(m4) 
Equating the terms up to fourth order in m one readily finds 

^2 ^ F2(A,0)-2dFo(A,0) 

rfF2(A,0)-|F4(A,0) • ^ ^ 

On Fig. [3] we plot the derivative of ln(m(A)^)/2 as computed from Eq. ffTSjl versus 
ln(A) with: 

A:=(A-Ac;)/Ac;. (19) 

In the limit A ^ we obtain (3 = 0.5 which is what should be expected for the 
universality class of the four dimensional Ising-model. 



4 Monte Carlo study 

In this section we study the model by means of Monte Carlo simulations. We first 
outline algorithms to sample field configurations and discuss their performance. In 
section 14.21 we present a phase diagram of the model. Subsequently, we analyze the 
mass gap of the theory at zero and finite temperature and estimate the relation 
between the bare coupling anisotropy 7 and the renormalized, physical anisotropy 
.^aniso ill sectiou 14.31 Finally, section 14.41 presents a finite-size scaling analysis where 
we determine Ac, several critical exponents, and the order of the phase transition. 
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4.1 Monte Carlo algorithms for the model 

Starting from Eqns. ([3]) and (jlj), the contribution to the overall action of a fixed site 
n can be written as 

S{^n) = -{Xn- ^^^Y + -{yn- ^^^^) ^ - W (x„ , . (20) 

That is, the action can be decomposed as 

5'[$] = S{^n) + (terms independent of$„) . (21) 

First, we implement a standard local Metropolis algorithm [S]. From Eq. (120]) one 
sees that the probability distribution for a single field variable in an otherwise 
fixed field configuration is given bj{§ 



p{^n) oc e 



-5(<I>„) 



OC exp{-^ (x„ - Xo)^} exp{-^ (y„ - yo)^} exp{ w(xn,yn)} • ^ ^ 

A local Metropolis algorithm can be implemented by generating proposals x' and 
y' for the new real and imaginary part of $„ from Gaussian distributions with the 
parameters found in eq. 0221) . The trial field variable = x' + iy' is accepted with 
probability min { 1 , exp[w(x', y') —w{xn, yn) ] } ■ We shall refer to this algorithm as 
M-G in the following. 

A different way to set up a local Metropolis sampler is to generate trial variables 
x' and y' from a homogeneous distribution on the intervals [xq — a^jXQ + a^] and 
[yo — ay,yo + ay], respectively. Here and ay are constants which can be chosen 
arbitrarily. Given the proposal <l>^ = x' +iy' one then accepts $' with the probability 
min{ 1 , exp [SE{^n) — Se{^')]}- We dub this procedure M-h. The efficiency of 
this algorithm will strongly depend on the choice of the interval widths ax and ay. 
It is expected that the algorithm M-h can be rendered more efficient by implementing 
it as 'modified Metropolis algorithm' [U] which performs several subsequent updates 
of the same site during each sweep over the lattice. We will test this procedure as 
well and refer to it as M-m. 

Since we anticipate the occurrence of a continuous phase transition (cf. section [3]) 
we also implement a Cluster algorithm in order to reduce the possible impact of crit- 
ical slowing down. We adapt the methods discussed in [TUl E] to our model. Our 
implementation of the cluster algorithm is outlined in detail in the appendix IA.2[ 
Given that the cluster algorithm only refreshes a subset of the degrees of freedom 
it has to be combined with one of the local Metropolis samplers described in the 
beginning of this section. We have chosen the algorithm M-G for this purpose and 
we refer to this combination of the cluster algorithm and M-G as GA in the following. 



^ w{x, y) has been defined in Eq. 
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Figure 4: (Color online) Plots of the order parameter and the susceptibility as 
obtained using the algorithms discussed in the text. The numbers in the legend are 
related to the algorithms in the following way: (1) M-G; (2) M-h; (3) M-m; (4) CA. 



We now report on our practical experiences with the above algorithms. In or- 
der to compare the algorithms we have carried out simulations using each of the 
algorithms on a 12^- lattice at fixed A = 0.5 varying the coupling anisotropy 
7. After initializing the lattice in a disordered way ('hot start') 20000 thermal- 
ization sweeps were performed and 1000 configurations were subsequently ana- 
lyzed, separated by 50 sweeps each. We have computed the expectation value of 
the field's modulu^ ( |$| ) := | J2n 1/^4 and the corresponding susceptibility 
X = V4 [( 1$ P ) — ( 1$ I )^ ] (V4 denotes the four dimensional volume of the lattice.). 
The results are displayed in Fig. HI As can be seen, there is a phase transition 
at 7 ~ 14 (This will be discussed in detail in section 14. 2[ ) Below this transition 
(7 < 12) the results obtained with all algorithms agree within the margins of errojfl. 
However, deviations exist for larger values of 7. We think that these are due to 
inefficiencies of some of the algorithms in the vicinity of the phase transition. The 
M-/iupdate performs particularly poorly in this range. The acceptance rate was very 
low and this could not be significantly improved by tuning the widths and ay of 
the proposal distribution. Its modified version M-m performs significantly better; 
indeed the data is in better agreement with the data generated by means of the 
other algorithms. We have decided to use the algorithms M-G and CA exclusively 
in the remainder of the work because they perform better and do not require tuning 
the widths of the proposal distributions. 

We close this section by briefiy describing some other technical details. In Tab.[T]an 
overview is given on the statistics and algorithms used to obtain the results quoted 
in this work. Here, rithermai denotes the number of thermalization sweeps that were 

^Cf. [12] for a discussion of the validity of this procedure and its systematic errors. 
^These were estimated by means of the jackknife method. 
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data displayed/quoted in 
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^'Sweeps 


'^config 
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Fig. \4\ 
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Sec. 112] 


10000 to 20000 
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1000 


M-G 


Tab. E] 
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M-G 


Tab. E] 
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M-G 


Tab. H 


10000 


50 


10000 (7 = 2), 20000 (7 = 3) 


M-G 


Sec. 113] 


10000 


50 
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CA 



Table 1: Summary of the statistics of the simulations in this work. (*): The number 
of configurations depends on the lattice size. - Confer the text for further explana- 
tions. 
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Figure 5: (Color online) Resolving the phase structure by running simulations cut- 
ting through the (7, A)- plane. This figure presents results obtained at fixed values 
of 7. 



carried out before observables were measured the first time after initializing the 
lattice. Tisweeps stands for the number of updates of the entire lattice carried out be- 
tween every two successive measurements, riconfig gives the number of configurations 
from which observables were calculated. We have tested ordered, disordered, and 
mixed phase initial configurations and checked that they yield the same results. We 
have always applied periodic boundary conditions. 

The margins of error quoted in this work and displayed in the figures were obtained 
by means of the jackknife method if not stated otherwise. 

4.2 Phase diagram of the theory 

In this section we present a phase diagram of the model in the plane spanned by 
the parameters A and 7. From the mean-field analysis (cf. section [3]) we expect 
spontaneous symmetry breaking to occur within a certain range of this plane. 
The phase structure is resolved by running simulations varying one of the param- 
eters while keeping the other one fixed. Fig. [S] shows results for the expectation 
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value of the field (i. e. the order parameter for the system under consideration) and 
the susceptibility as obtained from simulations where the coupling A was increased 
at different values of the anisotropy parameter 7. The sudden increase in each of 
the curves displayed in the upper panel together with the simultaneous peak in the 
susceptibility provides strong indication for a phase transition. Given that the re- 
sults for the order parameter $ appear to be continuous we suppose that the phase 
transitions are continuou^. From the peaks of the susceptibilities we have read off 
the critical values of the parameters and estimated their errors from the width of 
the peaks. Plotting the phase transition points in the (7, A) - plane yields Fig. [6] 
which provides the phase diagram of the model. Briefly summarizing its content in 
words, one can say that there is a line of presumably second-order phase transitions 
separating the symmetric phase at small A and/or high 7 from a symmetry broken 
phase which exists at larger values of A and not too high values of 7. Two features 
deserve to be commented on. First, it appears as if there is a critical Ac below 
which spontaneous symmetry breaking never occurs. This can be inferred from the 
results of the simulations keeping 7 = 1.0 fixecfl. Second, the numerical value of 
this critical coupling can be roughly estimated as Ac ~ 0.3 which agrees remarkably 
well with the prediction from the mean field analysis quoted in Eq. (ITSll I^. We notice 
that there are larger errors at high 7. This will be discussed in section [5l 

We briefly comment on the role of the parameter A = A ■ which has the di- 
mension of a mass. When investigating the phase structure of the lattice model we 
consider the dimensionless A as the relevant parameter. On the other hand, we can 
understand A as a cutoff that is externally given. By virtue of the relation between 
A and A we can then interpret as the spatial resolution or, equivalently, as a 
momentum cutoff. 



4.3 Mass gap and physical anisotropy 

In this section we determine the mass gap of the theory at both zero and finite 
temperature and analyze the dependence of the physical anisotropy .^aniso on the 
corresponding input parameter 7. 

First, we study the mass gap of the theory at zero temperature, i. e. on isotropic 
lattices where N^. = and 7 = 1. Before presenting our Monte Carlo results we 
provide an estimate for the mass from the second derivative of the potential at its 
minima. Keeping in mind that there is a factor 1/2 in front of V{^) in Eq. ([T]) the 

^We will provide evidence for the continuous character of the phase transition at 7 = 1 and 
7 = 10 in section [4.41 where we carry out a finite-size scaling analysis. 

^Physically, we interprcte the case 7 = 1 as the zero temperature limit while higher numerical 
values correspond to non-zero temperature. Thus we consider 1.0 as the lowest attainable value of 
7- 

^Also this issue is addressed in more detail in section l4l4l 
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Figure 6: Phase diagram of the model. 
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Table 2: Results for the masses at zero temperature, i. e. 7 = 1. The last column 
gives the tree level estimates in units of computed from Eq. ( l23i) . 



estimate for the mass at tree level is 



mo 



2 d{Re^y y 2 d{lm^y 



2A . 



(23) 



This will serve for later comparison to the numerical results. Notice that mg is as- 
sociated with the mass of topologically trivial fluctuations of the field $ about the 
minima = ±A of its potential. The fact that rriQ is larger than implies 
that these fluctuations do not take place. 

We have numerically determined the lowest mass of an intermediate state by mea- 
suring effective masses computed from timeslice correlators. It proved beneficial to 
run the calculations on lattices with different linear extensions because this allows 
to find a reasonable trade-off between systematic errors (due to too short a time 
axis) and statistical errors (due to limited computer time when simulating larger 
volumes). Tab. [2] presents the results for different choices of the coupling A. 
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Table 3: Results for the physical anisotropy. d in the fifth column denotes the 
separation at which the effective mass has been evaluated. 



In the next passage we present estimates for the physical anisotropy ^aniso = Os/flr 
resulting from a given bare anisotropic coupling 7. To this end, we proceed as follows. 
We require that the physical correlation length ^ be the same when measured along 
temporal and spatial axes. That is, 

e = a«6 = ttr^r (24) 

which implies 

Saniso ~7~ \^^) 

ttr is rris 

where we have used is = ^/'ms (and respectively for ir ) in the last equality, 
and rrir denote the masses measured along the respective axis in units of the inverse 
lattice spacing. 

We now make the assumption that .^aniso/7 = C'(l)- This was found to hold in the 
case of 0^- theory [13]. Under this assumption the lattice is approximately isotropic 
when the edge lengths are chosen as = 'jN^. This choice of the lattice extensions 
ensures that the impact of undesired finite temperature effects can be neglected. 
Using such a setup, we have measured effective masses as a function of distance and 
we have estimated ^aniso from 

Saniso — / 7\ • V^*-"/ 

rriri'jd) 

The purpose of evaluating rrir at separation 'yd is to ensure that the contributions of 
higher excitations have decayed away to the same extent when comparing the two 
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12 


12 
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Table 4: Results for the masses at finite temperature. 



effective masses. Tab. |3] summarizes our results from this procedure. As can be seen 
there, the physical anisotropy .^aniso agrees with 7 within the margin of error. In 
particular, this implies for the temperatur^ T = ^aniso/('2s iVr) = 7/(os A^r)- This 
shows that 7 is related to the physical temperature in a straightforward way and 
therefore it a posteriori justifies our use of 7 as a parameter in the phase diagram 
in section 14. 2[ It has to be noted though that we have only checked the relation 
between ^aniso and 7 for 7 < 3. We assume that we can extrapolate these results to 
larger values of 7. 

Finally, we have studied the temperature dependence of the mass. To this end, 
we have applied the same procedure as in the zero temperature case but this time 
running the simulation on an anisotropic lattice with 7 > 1. The results are dis- 
played in Tab. |H By comparing this result to the corresponding one in Tab. [2] one 
sees that the mass has slightly dropped but in general it does not exhibit a strong 
temperature dependence. We will comment on this in the conclusions (cf. section 

ED. 



4.4 Finite-size scaling analysis for the susceptibility 

In this section we carry out a finite-size scaling analysis of the Monte Carlo data for 
the susceptibility x- This will enable us to determine several critical exponents and 
hence the universality class of the model. Furthermore, we confirm that the phase 
transition when varying A at fixed temperature is indeed continuous. 
The data underlying the discussion in this section have been obtained from simula- 
tions on lattices of various linear extensions L. We have varied the coupling A at 
several choicef^ of 7aniso- In analogy to the definition of a reduced temperature we 
define a reduced coupling A as 

A :^ ^ (27) 

where Ac denotes the (a priori unknown) critical coupling. We denote the corre- 
lation length in the infinite volume limit by ^. Then, the standard finite-size scaling 



^We work in units where fc^ = 1. 

^"Deviating from our conventions in the rest of this work we will denote the bare anisotropy 
parameter by 7aniso in section [44] in order to avoid confusion with the critical exponent 7. 

^-^In fact, we should write this as Ac(7aniso) since we study the model at various, fixed choices 
of 7aniso- Howcvcr, wc will suppress the 7aniso- dependence in our notation for the sake of brevity. 
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niimljpr of data noints 

XX LXXXX V^X V_y X V_XCAj U i-Xi V-J XXX U 


Ar- 




ij 


/d o f 


1.0 


51 


0.3299 (3) 


0.98 (3) 


0.50 (2) 


2.2 


2.0 


140 


0.3719 (3) 


0.924 (13) 


0.449 (11) 


1.7 


5.0 


124 


0.3810 (7) 


0.84 (3) 


0.421 (14) 


1.6 


10.0 


103 


0.352 (7) 


1.05 (6) 


0.67 (5) 


1.2 



Table 5: Results for Ac and the critical exponents 7 and v at several values of 7amso 
as obtained from the finite-size scaling method. The last column gives the value of 
per degree of freedom. The margins of error are the estimates as obtained from 
the standard x^- fit- Confer the text for a comment on systematic errors. 



approach states that the susceptibility Xhi}^ as measured on a finite box of linear 
extension L is given by 

XL(A)=r/'^Xo(Ve). (28) 

Here, 7 and v denote the critical exponents such that in the thermodynamic limit 
X ~ A~'^ and ^ ~ in the vicinity of the critical point at A = 0. Xo is a continuous 
function which only depends on the ratio L/^. Substituting the asymptotic behavior 
of ^ into the last equation and making the ansatz 

X,{Ly)^{^^''^ !{L'I^\) (29) 

with some continuous function / one arrives at 

Xl(A)=L^/V(^'/^A). (30) 

Assuming an expression for / one can fit data from lattices of different lengths using 
the last equation and thus one can determine the critical exponents and A^ (The 
dependence of Eq. ( l30l) on Ac enters via the definition of A in Eq. (|27j) .). 
Introducing a set of six fit parameters ao to our ansatz for fitting / reads 

y(k,L) ■= L"^/"^ . (31) 

^l + al[LiM(A_i)_a^]2 

The significance of the fit parameters is as follows, ao corresponds to the value of the 
critical coupling Ac, a\ and 02 correspond to the critical exponents 7 and respec- 
tively, as serves to take into account that the peak of the susceptibility as measured 
on a finite volume occurs slightly shifted away from its position at A = in the 
thermodynamic limit. The parameters 03 and 04 are not of any immediate physical 
significance but serve to parameterize the function which resembles a Lorentz curve. 
Tab. [5] summarizes the results from fitting our Monte Carlo data while Fig. [7] visu- 
alizes the outcome of the finite-size scaling procedure for 7 = 1 and 7 = 10. Several 
things need to be commented on. First, the margins of error quoted in Tab. EJare 
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Figure 7: (Color online) Test of the finite-size scaling. The values inserted for Kci 
7, and z/ are the ones quoted in Tab. O 



the statistical errors obtained from the x^- fito in the common way; they appear to 
be underestimated, in particular for Ac. We think that there are systematic errors 
that have not been taken into account in Tab. E] which affect the fit. This is be- 
cause the statistical errors of the Monte Carlo data are largest in the vicinity of the 
critical point where the finite-size scaling method should make the curves collapse. 
Therefore, the x^- fit gives a lower weight to the data that are more relevant for 
physical reasons. It is difficult to estimate the numerical size of this effect. 
Comparing the quality of the fit at the different values of 7aniso (cf. Fig. [7]), one sees 
that it works better at smaller 7aniso- In the limit of high temperature (Fig. [Tl^b) ) 



^^X^ is the common merit function of the x^- fit and not the square of the susceptibility x- Since 
there should be no risk of confusion we do not introduce an extra piece of notation to distinguish 
the two quantities. 
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the statistical errors are significantly larger. In particular, we think that one should 
be a bit careful about the estimate for Ac(7aniso = 10) ~ 0.35 which appears unre- 
liable when compared to the value read off from Fig. [5l 

Giving an interpretation of these results, we first note that the critical exponents 
7 and v obtained at 7amso ^ 2 agree well with the ones for the Ising model in four 
dimensions. Thus we deduce that the model lies in the same universality class as the 
4D-Ising model for these choices of the parameters. At this point it is necessary to 
discuss the non-renormalizability of the potential ([2]). From an expansion about one 
of its minimal^ one finds that the potential contains contributions up to arbitrary 
order in the real and imaginary part of $. The corresponding coupling constant for 
a contribution of order n is proportional to A"'^""^-'. The terms with a coupling of 
negative mass dimension are irrelevant for the critical properties of the model [1] . 
The only relevant and marginal contributions are the terms of order and <I>^ 

and the critical exponents should be completely determined by them. 
The data at higher temperatures is not interpreted in as straightforward a way. 
From Tab. EJone sees that the critical exponent v shows a noticeable increase when 
7aniso is raised to 7amso = 10 after dropping at intermediate values. In the limit of 
high temperature (i. e. 7aniso = 10), the numerical value of v is compatible with the 
corresponding critical exponent of the three dimensional Ising model whose value 
is 0.630 ± 0.0015 [13]. This observation might hint at the dimensional reduction of 
the four dimensional field theory to a three dimensional one in the limit of infinite 
temperature (or equivalently 7amso — ^ cxo). Indeed, it was found in [151 [13 [13 [E] 
that a dimensional crossover can be observed in the critical properties of a field 
theory when the temperature is increased from T = towards T — oo. 

From the values quoted in Tab. we also infer that the phase transition is always 
a continuous one. For a first-order phase transition the susceptibility at the critical 
point diverges proportional to the volume [IHl EDI |2T] . From Eq. (l30i) it follows that 
X ~ {L/^)'^/'^. Given that our data is certainly incompatible with 7/// = 3 for all 
values of 7amso under consideration we can exclude that the phase transition is of 
first order. 

It is also interesting to numerically determine the exponent (3 which is defined by 
Eq. (fT6l) and which has previously been estimated in the mean-field approximation 
(see Fig. [2]). The result of a linear fit to ln($) plotted versus In A at 7aniso = 1 
is shown in Fig. [HI We have only fitted the data from the 24^- lattice because 
it provides the best realization of the thermodynamic limit among the available 
simulation results. The fit yields 

(3 = 0.453 ± 0.002 (32) 
where the error estimate gives the statistical error as computed from a x^- fit again. 

^•^Due to the pole at $ = such an expansion will only be valid in a sufRciently small vicinity 
of the respective minimum. 
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Figure 8: Double logarithmic plot of ($) vs. A. 



This result deviates about ten per cent from the MFA- estimate quoted in sect. [3l 
We emphasize that this way to determine (3 is prone to systematic errors which have 
not been taken into account in the above error estimate. Moreover, we expect that 
our approach to fit /3 systematically produces too small a value because the curve 
of ($) will always be a smooth one when computed on a finite lattice. 

We conclude this section by presenting the Binder cumulant for the zero temperature 
case (i. e. 7aniso = 1) in Fig. [9l The Binder cumulant [22] is defined as 

1 

'''' 

In Fig. Pthe curves intersect at Ac — 0.33, which is in reasonable agreement both 
with the finite-size scaling result quoted in Tab. [5] and the mean-field estimate in 
Eq. ( fT5|) . Furthermore, the numerical value of the Binder cumulant at the critical 
point is characteristic for each universality class. In [23] it was found by analytical 
means that for the four dimensional Ising model (s^)/ (s^)^ = 2.188... where s is the 
average magnetization per unit volume. Numerical studies of the four dimensional 
Ising model [2^ [25] have yielded = 1.92(3). The corresponding numerical 

values for the Binder cumulant U{L,Ac) as defined in Eq. (l33l) are 0.27 and 0.36 
(independent of L, of course), respectively, and they are indicated in Fig. [9] by the 
solid and dashed line. It seems as if our data favors {s^) / {s"^)^ = 1.92. However, 
the data does not allow for a completely unambiguous discrimination against one of 
these two alternatives. Nevertheless, we infer that the findings for U{L, Ac) support 
the statement that the model at hand at zero temperature falls into the universality 
class of the 4-d Ising model. 
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Figure 9: (Color online) The Binder cumulant U{L,A) plotted vs. A as computed 
on lattices of various sizes. The lines connecting the data merely serve to guide the 
eye. The dashed horizontal line corresponds to the result quoted in [23] while the 
solid line corresponds to the results of [211 EH] • 



5 Summary and conclusions 

In this work we have analyzed the theory governed by the Lagrangian ([T]) by both 
mean-field and Monte Carlo methods. At zero temperature, the mean-field analysis 
predicts the existence of a continuous phase transition and the universality class 
of the four-dimensional Ising model. Both results are in agreement with the Monte 
Carlo results. It seems that the additional symmetry of the potential V{^) = V{^*) 
does not have any impact on the critical behavior of the system. 
Apart from the phase structure, we have also determined the mass gap of the theory 
both at zero and finite temperature. We found that its value is almost independent 
of temperature. The physical anisotropy did not show deviations from the bare 
coupling anisotropy 7 within the margin of error. Finally, we obtained the critical 
exponents for the transition when varying A also at non-zero temperature from a 
finite-size scaling analysis. In the high temperature limit, the numerical values seem 
to indicate dimensional reduction to the three dimensional Ising model. However, 
more work is necessary to understand the situation at high temperature. 
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A Appendix 



A.l Some formulae for the MFA 

In this appendix we summarize some technical details of the mean-field approxi- 
mation presented in section [SI We use 7 = 1. It is straightforward to extend the 
analysis to 7 > 1, however, only for 7 less than about 2 are the results, at least 
qualitatively, reasonable, for larger anisotropy a more involved analysis would be 
necessary. 

We start by noting that the functions F„(A, m) defined in Eq. (11 2p admit the fol- 
lowing series expansion in m: 

FniKm) = {2dY^-{2dmf e-^^^M rfe(cose)"+'=e^'^°^(2®) 

fc=o ■ -^^ 

(34) 

Note that the terms with odd n + k vanish by virtue of the symmetry properties of 
the 0- integral. From this expression we can obtain the derivatives at m = 0, to be 
used in ([HI), f[T7|l . in terms of the integrals: 



i?„(A) = - / di?i?2n+ig-v(i?) ^ 2(^^) ' Kn+i{A'\^2d + k^) (35) 

= i f ' de(co==e)"e*'-w = ^ |; ( ^ ) ^i„(„)|„^,,(36) 

The uj = ^ solution can be obtained by reversing the sign of the cos(20) term in 
([MD. 

A. 2 Implementation of the cluster algorithm 

This appendix describes our implementation of a cluster algorithm for the model. 
We start by rewriting the lattice action Eq. ([3]) in terms of the real and imaginary 
parts Xn and y„ of each field variable $„, which yields 

n j 

+ 2(3 + 72)11)^2^^(1.,)}. 

We now identify a set of Ising- like degrees of freedom in the model by defining 
Sn '■= sign(x„) = sign (Re$„) which we will use to set up the cluster algorithm. 
The action can then be rewritten as S" = S'ising + S'Rest where 5'Rcst does not depend 
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on the {sn}n and 
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^ ^ ^ ^ Jn,n+fi Sn Sn+fj, 
n fj, 



(37) 



__1 V- , 

<l,n> 

Here we have defined a set of site-dependent nearest neighbor couphngs as 

{7 \xr,\ if 1, n are time-hke neighbors 
iL L • (38) 

- \xi\ \xn\ if 1, n are space-hke neighbors 

We have chosen to embed Ising spins only in the real part of the field. Due to the 
shape of the potential we expect the imaginary part to be less affected by critical 
slowing down, if at all. Moreover, we note that it is not possible to define Ising 
spins from a projection on a random axis in the complex $„- plane as it can be done 
for instance in the case of the XY- model. The reason for this is that the potential 
does not possess the necessary symmetry; it is in general not possible to fiip the 
projection of $„ on an arbitrary axis without altering the value of S'Rest- 
It remains to set up a cluster algorithm for the Ising- like degrees of freedom. The 
dynamics of the spin variables is equivalent to the one of a four dimensional 

Ising model governed by the Hamiltonian 

Hlsmg = 'S'ising = — ^ ^ Jl^n ^1 ■ (39) 
<l,n> 

We then proceed in the standard way [291 [3Ql [31] by introducing bond variables di^n 
for all pairs of nearest neighbors /, n. From the literature (cf. e. g. [30]) it is known 
that the probability for a bond to be active has to be chosen as 

- exp(-2J/^„) if si = Sn 

in order to implement the correct dynamics for the Ising Hamiltonian ( l39i) . After 
activating each bond variable with this probability we proceed as described in [31] : 
First, we choose a lattice site at random. Then we track the percolation cluster to 
which this site belongs and fiip the real parts of all field variables in the 

percolation cluster. 
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